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A study of the formation of excitons as a problem of two Dirac particles in a gapped graphene 
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limit the separation of the center-of-mass and relative motions is obtained. Analytical solutions 
for the wave function and energy dispersion for both cases when electron and hole interact via a 
Coulomb potential are found. It is shown that the energy spectrum and the effective exciton mass 
are functions of the energy gaps as well as interlayer separation in case of two layer gapped graphene. 
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I. INTRODUCTION 



The exciton formed by an electron and a hole is one of the very important objects in both experimen- 
tal and theoretical physics of semiconductor heterostructures and graphene. The study of the excitonic 
system and its properties is necessary because of potential applications in electronics and photonics, 
including, design of thresholdlcss lasers, optical computing and quantum computing [l|-IH- From the 
theoretical point of view, an exciton is a two-body system and to address the problem we have to solve 
a two-body problem in semiconductor heterostructures or graphene. Let us mention that the problem of 
the interaction between two particles is very important for the deep analysis of the many-body physics 
in excitonic system. The analysis of the spectrum of collective excitations necessary to study such many- 
body phenomena in the excitonic system as Bose-Einstein condensation and superfluidity 0] requires 
the solution of two-body problem for a single exciton. The temperature of Kosterlitz-Thouless phase 
transition corresponding to the superfluidity depends on the density of the superfluid component (Io| . 
This density of the superfluid component is defined by the two-particle Green function of weakly inter- 
acting excitons. The problem to find the two-particle Green function of the interacting excitons requires 
one-particle Green function of non-interacting excitons, which is defined by the wave function and energy 



dispersion of a single exciton [11[ . The collective properties and superfluidity of excitons and polaritons 



tive prop 

in gapped graphene have been studied in Refs.fl2l.ll3l Therefore, it is of fundamental and practical inter- 
est to focus on solution of two-body problem in semiconductor heterostructures and graphene. Finding 
solution of two-body problem in gapped graphene layers is the subject of this Paper. 

To describe the formation of an exciton in semiconductor heterostructures like quantum wells the 
standard quantum mechanical approach is used based on Schrddinger equation. In this case the two- 
body problem with a scalar action-at-a-distance inter-particle potentials is completely understood and 
dev elop ed both in configuration or momentum three dimensional as well as two dimensional (2D) spaces 
[Tj-Qjl . The notions of absolute time and absolute space allow us to describe the two particles of masses 
mi and m 2 with Euclidean position vectors ri and r 2 , and momenta pi and P2 in an inertial frame. It is 
well known [lij that in the covariant mechanics, the most straightforward solution of the classical problem 
with scalar interaction is obtained by a transformation from the individual particle coordinates ri , r 2 
and pi , p2 to the covariant center-of-mass and relative coordinates. Using canonical transformation 



m, 1 r 1 +m 2 r 2 

R = , r = n r 2 and P = pi + p 2 , p = pi p 2 , (1) 

m 

where m = mi + m 2 is a total mass of the system, we can separate the decoupled center-of-mass from 
the relative motion in configuration or momentum space. In these coordinates, the nontrivial motion of 
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the system of two particles occurs entirely in the reduced problem of one body motion and motion of the 
center-of-mass. The total Hamiltonian of the system can be presented as a sum of two parts. One part 

p2 

describes the motion of the center-of-mass while the other two particles relative motion: H=-^+ H r 

2 

with H r = + V(r), where [i = mimz/m is two body's reduced mass. The Hamiltonian H r governs 
the relative motion of two particles and, when the Schrodinger equation with this Hamiltonian have been 
solved, the wave function of the relative motion of two particles is obtained. 

Today graphene has been attracting a great deal of experimental and theoretical attention because of 
unusual properties in its band structure 17. 18|. Graphene is a two-dimensional layer of carbon atoms, 
where the atoms form honeycomb lattice (l9l l2dl|. In the low-energy limit the low energy excitations 
in graphene are described not by the Schrodinger equation but instead in graphene electron and holes 
behave as relativistic massless particles described by a Dirac-like equation for masslcss and chiral particles 
[2l| - |25j , which is known as Weyl equation. Many of the unusual properties of graphene arise from the fact 
that its quasiparticles are described by Dirac spinors. Since electron and holes in graphene are governed 
by the Weyl equation they have an intrinsic degree of freedom that resembles the spin degree of freedom 
in the original Weyl equation. This degree of freedom is called pseudospin in order to distinguish it 
from the spin and is described by the Pauli matrices. In connection to the pseudospin, there is a good 
quantum number - the chirality that is defined to be the projection of the 2-momcntum operator on the 
direction of the pseudospin. Clearly, the electrons will have positive chirality and the holes will have 
negative chirality. Also the Fermi speed vf ~ 10 6 m/s that is around 300 times slower than speed of 
light replaces the speed of light in the original Weyl equation. The positive-energy solutions of this 2D 
equation describes electrons, whereas the negative-energy solutions describes holes. 

A 2D exciton which is a bound state of an interacting electron and hole in gapped graphene presents a 
two-body system and it is a fundamental and practical interest to focus on solution of two-body problem 
in graphene. Graphene consists of two equivalent carbon sublatticcs and quantum-mechanical hopping 
between the sublatticcs leads to the formation of two energy bands, and their intersection near the edges 
of the Brillouin zone yields the conical energy spectrum. As a result, quasiparticles in graphene exhibit a 
linear dispersion relation. The formation of excitons requires the existence of the energy gap in electron 
and hole energy spectrum. The formation of excitons in gapped graphene was studied in Refs. [l2l [26l - [28j . 
The electronic ground state of intrinsic graphene and bilayer graphene in the absence of the energy gap 
using density functional theory within the local-density approximation and Bcthe-Salpeter equation was 
studied in Ref. [2^] . According to Ref. [2!| , no pure bound exciton was identified in intrinsic graphene and 
bilayer graphene in the absence of the energy gap. Therefore, excitons in graphene can be formed due to an 
energy gap opening in the electron and hole spectra in the graphene layer. There are different mechanisms 
of electronic excitations in graphene. The energy gap in graphene can be induced and controlled by a 
magnetic field, doping, an electric field in biased graphene, and hydrogenation [26L [30l - 432j ] . 

By contrast quantum mechanics in the description of relativistic two-body problem the situation is 
extremely more complicated and till now there is no completely self-consistent theory for the separation 
of the center-of-mass and relative motions even for the two-body case. This is due to the following facts: 

i) the particles locations and momenta are 4-vectors, 

ii) the momenta are not independent but must satisfy mass-shell conditions, 

iii) the inter-particle interaction potentials appear in the boosts as well as in the energy generator 
in the instant form of dynamics. As a result of a transformation to the center-of-mass system even a 
scalar action-at-a-distance inter-particle interaction potential becomes dependent on a coordinate and 
momentum. 

iv) the structure of the Poincare' group implies that there is no definition of relativistic 4-center of 
mass sharing all the properties of the non- relativistic 3-center of mass [33| . 

The two body problem in graphene when the electron and hole are interacted via a scalar action-at- 
a-distance inter-particle potential and are described by Weyl's equation becomes even more complicated 
than in simple relativistic case. Firstly, it is related to the fact that speed of light in the Weyl equation 
for graphene is replaced by the the Fermi speed and the resultant equation becomes non-covariant and 
canonical transformation implementing the separation of the center-of-mass from the relative variables 
within relativistic approach is invalided. Secondly, even the inter-particle interaction depends only on 
the coordinate of the relative motion after to the center-of-mass transformation it becomes dependent 
also on a momentum. Lastly, the center-of-mass energy can not be separated from the relative motion 
even though the interaction depends only on the coordinate of the relative motion. This is caused by the 
chiral nature of Dirac electron in graphene. 



3 



Our paper is organized in the following way. In Sec. |TT] we present the Hamiltonian of the spatially 
separated electron and hole in two different parallel graphene layers separated by a dielectric in the 
presence of the band gap. In Sec. Illll is given the procedure of separation of the center-of-mass and relative 
motions for two particles in graphene. In Sec. IIVI we obtain the energy spectrum and wave function of 
dipole exciton in two-layers graphene and find the effective exciton mass. The energy spectrum and wave 
function of exciton formed in a gapped graphene layer and corresponding effective exciton mass arc given 
in Sec. |V] Finally, the conclusions follow in Sec. IVI1 



II. EXCITON HAMILTONIAN 



Let us consider two different parallel graphene layers with the interlayer separation D and assume that 
excitons in this system are formed by the electrons located in the one graphene layer and, correspondingly, 
the holes located in the other. In this system electrons and holes move in two separate layers with 
honeycomb lattice structure. Since the motion of the electron is restricted in one graphene layer and 
the motion of the hole is restricted in the other graphene layer, we replace the coordinate vectors of the 
electron and hole by their projections ri and r2 on plane of one of the graphene sheet. These new in- 
plane coordinates ri and T2 will be used everywhere below. Thus, we reduced the restricted 3D two-body 
problem to the 2D two-body problem on the graphene plane. Each honeycomb lattice is characterized by 
the coordinates (tj, 1) on sublattice A and (rj, 2) on sublattice B with j = 1,2 referring to the two sheets. 
Then the two-particle wave function, describing two particles in different sheets, reads ^(ri, si; r 2 , s 2 ), 
where ri and r 2 represent the coordinates of the electron and hole, correspondingly, and si, s 2 are 
sublattice indices. This wave function can also be understood as a four-component spinor, where the 
spinor components refer to the four possible values of the sublattice indicies si, s 2 : 



*(rx,si;r 2 ,s 2 ) 



/ 



V 



b aa (ri,r 2 ) 
*>a&(ri,r 2 ) 
*>6a(ri,r 2 ) 
*Wri,r 2 ) 



where ^> Q = 



'a a 
t>ab 



(2) 



The two components mean that one particle is on sublattice a(b) and the other particle is on sublattice 
a(b), correspondingly. In other words, the spinor components are from the same tight-binding wave 
function at different sites. Each graphene layer has an energy gap. Obviously the energy gaps in graphene 
layers are independent and in the general case we can introduce two non-equal gaps S± and 62 for the 
first and the second graphene layers, respectively. 

The corresponding hopping matrix for two non- interacting particles, including the energy gaps S\ and 
62 on the first and second layers, correspondingly, then reads 



Si 



Hn — 



V 



4 
4 




d 2 di 

-S 1 - 6 2 di 
£1+62 d 2 



/t 



zt 



(3) 



S 1 -S 2 J 



([3]) g?i = Kvp(—id Xl — d yi ), c?2 = JiVF(—id X2 — d V2 ) and the corresponding hermitian conjugates 
hvp{—id Xl + d yi ), 4 = h v F{—id X2 + d V2 ), where d x = d/dx and d y = d/dy, x\, y\ and X2, 
2/2 are the coordinates of vectors ri and r2, correspondingly, and vf is the Fermi velocity of electrons in 
graphene. This Hamiltonian allows us to write the eigenvalue equation for two non-interacting particles 



In Eq. 

are d\ = 



H & = (0 ^0 



which leads to the following eigenenergies: 



(4) 



e (Jfci, 5i; k 2 , S 2 ) = ±\jh 2 kl + m\ + <5? + 8\ ± 2^ \tfk\ + 6f)(h 2 k% + 5 2 2 ) = ±^ h 2 k\ + S\±^r9k\ + Sj . 

(5) 

where k\ and k2 are momentum of each particle, correspondingly. Eq. ([5]) gives the energy spectrum 
for two non- interacting particles in the presents of the non-equal gaps energies Si and 62 ■ The energy 
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dispersion is symmetrical with respect to the replacement of particles 1 and 2. When there are no gaps, 
Si = and S 2 = 0, as it follows from §5§ the energy dispersion is ±hk\ ± Kk 2 . 

Let's consider the electron and hole located in two graphene sheets with the interlayer separation D, 
and interacting via the Coulomb potential V(r), where r is the projection of the distance between an 
electron and a hole on the plane parallel to the graphene layers. Now the problem for the two interacting 
particles located in different graphene layers with the broken sublattice symmetry in each layer can be 
described by the Hamiltonian 

/ -6 1 + S 2 + V{r) d 2 dx \ 

4 S!-S 2 + V(r) d t 

d\ 8 1 +S 2 + V(r) d 2 

V d\ 4 S 1 -S 2 + V(r)J 

and the eigenvalue problem for Hamiltonian (|6|) is 

= e#, (7) 

where 'J are four-component eigenfunctions as given in Eq.@. 

The Hamiltonian (J6]) describes two interacting particles located in two graphene layers and satisfies the 
following conditions: 

i) when both gaps vanish Si = and 5 2 = 0, as well as two-body potential V(r) = 0, the Hamilto- 
nian describes two non-interacting Dirac particles. It is important to mentioned that eigenenergies are 
symmetrical with respect to the replacement of particle 1 and 2. 

ii) when the interaction between particles vanishes, V(r) = 0, it describes two independent particles, 
each located in the separate graphene layer, having two independent gaps energies related to the broken 
sublattice symmetry in each graphene sheet; 

iii) when the gaps in each graphene layer vanish, <$i = and S 2 — 0, the Hamiltonian describes two 
interacting particles in one graphene layer interacting via Coulomb potential V(r) and is identical to the 
Hamiltonian (2) in Ref. [27j representing the two-body problem in one graphene layer if the band gap is 
absent; 

iv) when the gaps Si = S 2 = S, the Hamiltonian describes two interacting particles in one graphene 
layer interacting via Coulomb potential V(r) = e 2 /4:ireoer, where e is the charge of the electron and £ is 
the dielectric constant of the graphene layer. 

v) when an electron and hole located in two different graphene sheets with the interlayer separation D, 
they interacting via the potential V(r) = — e 2 /AirsoedVr 2 + -D 2 , where £<f is the dielectric constant of the 
dielectric between two graphene layers. Let us mention that for 5i = S 2 = and D = the Hamiltonian 
([6]) is also identical to the Hamiltonian (2) in Ref. (27j . 



(6) 



III. SEPARATION OF THE CENTER-OF-MASS AND RELATIVE MOTIONS 

In Hamiltonian ([6]) the center-of-mass energy can not be separated from the relative motion even 
though the interaction V = V(r) depends only on the coordinate of the relative motion. This is caused 
by the chiral nature of Dirac electron in graphene. The similar conclusion was made for the two-particle 
problem in graphene in Ref. [27j . 

Since the electron-hole Coulomb interaction depends only on the relative coordinate, we introduce the 
new "center-of-mass" coordinates in the plane of a graphene sheet: 

R = avi + fir 2 , 
r = n - r 2 . (8) 

Here the coefficients a and /3 are to be determined later. Apparently we can use the analogy of the two- 
particle problem for Dirac particles in gapped two-layer graphene with the center-of-mass coordinates 
for the case of Schrodinger equation. The coefficients a and (3 will be found below from the condition 
of the separation of the coordinates of the center-of-mass and relative motion in the Hamiltonian in the 
one-dimensional "scalar" equation determining the corresponding component of the wave function. 
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To find the solution of ([7]) we make Anzatze 

^•(R,r) = e iK ' R ^(r) ■ (9) 
Let us introduce the following notations: 

K.+ = /Ce* 9 = K-x + ilCy , 
/C_ = /Ce 10 = K x — ilCy , 
-l I 



e = tan- 1 I ^ 1 , (10) 

and rewrite the Hamiltonian © in terms of the representation of the coordinates R and r in a form of 
the 2x2 matrix as 

v= ( 2 + V{r)a - 5ia + S 2 a 3 0i \ , . 

\ C\ 2 + V(r)a - 5kt + fcoa ) ' 1 > 

where 0\ and 02 are given by 

1 = hv F {aK--id x -d v )a Q , (12) 

«-*'(^A-*^ + o* + *) M 

where a; and y are the components of vector r, Uj are the Pauli matrices, o"o is the 2x2 unit matrix. 
Analysis of the operators (fT2|) and |13]) shows that the coordinates of the center-of-mass and relative 
motion can be separated in a certain approximation. 

For the Hamiltonian (|TTj) the eigenvalue problem = e^> results in the following equations: 

(02 + V '(r)<70- <Vo + <S 2 cr 3 )* a + 01*6 = ecr *a 

Cl* Q + (02 + V(r)<7 - *iod + S 2 a 3 ) = eo- * 6 • (14) 

From Eq. (jT^J) we have: 

* h = (ea - 02 - V(r)a + S^o - 8 2 a z y 1 0{V a . (15) 

Assuming the interaction potential and both relative and center-of-mass kinetic energies are small com- 
pared to the gaps 5± and S 2 we use the following approximation: 

(eer - 02 - V(r)<j + o"i<7 - S^y 1 ^ — — . (16) 

eao + dido — o 2 <73 

Using the fact that the operator 0\0\ is purely hermitian, applying Eq. (| 14[) and 

0{Oi = h 2 v 2 F (a 2 JC 2 - V 2 r - 2ia{K x d y + K y d x )) a , (17) 

we obtain: 

2 + V{r)a - <Vo + S 2 a 3 ) <J a + h 2 v 2 F ± r V * y -^-* a = ea V a . (18) 

eoo + (>i(Tq - d 2 a 3 

Now using Eq. ([2]) we rewrite Eq. (|18j) for the individual spinor components in the following form: 

- Sl + s 2+ V (r) + fiV /^-y;-y^+^ 

e — oi — 02 

?iu_F (/?/C_ + i5 x + ab = e0 oo , (19) 

hv F ( /3/C+ + id x - d y ) 4> aa + 

-S 1 -S 2 + V(r) + h 2 v 2 F £ f * 2-^ ab = e^ ah . 20 

e - o 1 + o 2 I 



G 



We solve Eq. ([20]) with respect to <j) ab : 
4>ab= e + 5i+ S 2 - V{r) - h 2 v F 



2 a 2 JC 2 -V 2 - 2ia(tC x d x + fC y d y ) 



e - Si + S 2 

Substituting </> Q b from Eq. (|2T|) into Eq. (|T9| . we obtain: 



( /3/C+ + id x - d y ) hv F (f> aa ■ (21) 



-<5i + 5 2 + V(r) + h 2 v 



2 2 o?K? -V 2 - 2ia(fC x d x + JCydy) 



6 — 81 — 82 



+ h 2 v 2 F {ptC- + id x + d y ) 
x ( /3/C+ + id x - d y ) = eq 



e + S 1 + S 2 - V(r) - HV F a !£ Z Yl Z «g ^ + *W 

e — 01 + 02 



(22) 



Assuming again that the interaction potential and both relative and center-of-mass kinetic energies are 
small compared to the gaps 8% and S 2 we apply to Eq. (1221) the following approximation: 



e + 5i+5 2 - V(r) - h 2 v\ 



a 2 fC 2 - - 2ia{tC x d x + JC v d v ) 



e — 8% + S 2 



1 



e + Si + S 2 



(23) 



Applying the approximation given by Eq. (|23|) to Eq. (f22|) . we get from Eq. (|22|) the eigenvalue equation 
spinor component 4> aa in the form: 



-ff! + S 2 + V(r) + ^ ^--V^-^A+M,) 



, 2 2 W - V? - 2i/3(/C X 9 X + JCydy) 

~ n V F J x 

e — 01 — S 2 



e — Si — 5 2 

aa = £<i 



(24) 



Choosing the values for the coefficients a and j3 to separate the coordinates of the center-of-mass (the 
wave vector JC) and relative motion r in Eq. (|24p . we have 



P = 



e — 5i — 5 2 
2e 

e + + J 2 
2e 



(25) 



After substitution of Eq. (f25j) into Eq. ([24|) we can obtain the component <j> aa of the spinor J2} as a 
solution of one-dimensional second order differential equation: 



{hv F JCf 



V(r) 



e(hv F ) 2 V 2 



<j>aa = [e+Sl — S 2 ] <j) aa 



2e ' ' v/ 2( £ 2 _ ( 5l +5 2 )2) 
The other components of ([2]) can be found as: 

* b = -(e<7 + (d X2 ai - idy 2 a 2 ) - Si<r - S 2 a 3 - ^(r)<T ) _1 iD\v a 
with Pauli matrices aj and 2x2 unit matrix <tq. Moreover, we have 

1 



e + Si+S 2 - V(r) + 



e-Si + S; 



(id X2 - d y 



(26) 



(27) 



(28) 
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IV. TWO BODY PROBLEM IN TWO GAPPED GRAPHENE LAYERS 



Let us consider an electron and hole located in two different parallel graphenc layers and interacting 
via the potential V(r) = — e 2 /AweoEdVr 2 + D 2 . Substituting this potential into (|26|) we obtain the second 
order differential equation for the component 4> aa . This equation cannot be solved analytically. However, 
assuming r <C D we can approximate V(r) by the first two terms of Taylor series and substituting 
V(r) = —Vq + 7r 2 , where Vq = e 2 / 'Airs^edD and 7 = e 2 /87reo£d-D 3 for the interaction potential into 
Eq. PS]), we obtain 



2e(hv F yVi 2 



(«i + s 2 y 



+ S 1 -S 2 + V a 



(hvplC) 

2t~ 



(29) 



The last equation can be rewritten in the form of the two-dimensional isotropic harmonic oscillator: 

(- J"i(e)V 2 + 7 r 2 ) <t> aa = F {e)4>aa , (30) 

where 

2e(hv F ) 2 



c 



2 - (Si + S 2 ) 2 ' 



T = e + 5l -S 2 + V ~ { -^. (31) 
The eigenfunction and eigenenergy for two-dimensional isotropic parabolic well were first determined 



by Fock in 1928 [341 ], later in Ref. [35[ , and was studied in detail in Ref. |36j. The single-particle 
eigenfunction for the two-dimensional oscillator widely used for the description of quantum dot [3 71 ] . The 
solution of Eq. (f3T)]) is well known (see, for example, Ref. [38]) and is given by 

2Ni rr 



•Fi(e) V 7i(e) ' 

where N = 2N + \L\ + 1, and N = min(n, n 1 ), L = n — ri', ri, n' = 0, 1, 2, 3, . . . are the quantum numbers 
of the 2D harmonic oscillator. The corresponding 2D wave function in terms of associated Laguerre 
polynomials can be written as 

«Wr) - ^^b^L)^-^^ x Lfir 2 ,^))^ , (32) 

1/2 




where 4> is the polar angle, L p k (x) arc the associated Laguerre polynomials and a = y\/F\{e)/ (2^/7) 

After some straightforward but lengthy calculations and the expansion up to second order in /C we 
obtain the following expression for the energy in quadratic order with respect to K, 

{hv F IC) 2 , (33) 

where fj, = 61 + S 2 and C\ = 2^N 2 (Hvf) 2 ■ Thus, from (|33[) we can conclude that the effective exciton 
mass M is given as a function of total energy gap S\ + S 2 as 

(34) 

The dependence of the effective exciton mass M defined by Eq. (|34| on the total energy gap <5i + 5 2 and 
the interlayer separation D for two graphene layers separated by the dielectric GaAs is shown in Fig. [TJ 
Firstly, according to Fig. [TJ the effective exciton mass M increases when the total energy gap Si + S 2 
and the interlayer separation D increase. Secondly, the mass increases much faster for the bigger value 
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of the interlayer separation and bigger value of the energy gap. Let's mentioned that the dependence 
of the effective exciton mass M on the interlayer separation D is caused by the quasi-relativistic Dirac 
Hamiltonian of the gapped electrons and holes in graphene layers. However, for the excitons in couple 
quantum wells (CQW) the effective exciton mass does not depend on the interlayer separation, because 
the electrons and holes in CQW's are described by a Schrodinger Hamiltonian, while excitons in two 
graphene layers are described by the Dirac- like Hamiltonian ((_>]). 




Value of Si+S 2 
— 100 meV 
-- 150 meV 
— 250 meV 
350 meV 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.8 0.9 1 1.1 1.2 1.3 1.4 

(5j+__ (eV) Interlayer separation D (nm) 

(a) (b) 



FIG. 1: Excitons in two graphene layers separated by the dielectric GaAs. The effective exciton mass M in 
the units of free electron mass mo (a) as a function of the total energy gap for the different graphene interlayer 
separations, (b) as the function of on interlayer separation for different values of the total energy gap. 



In Fig. [__ is shown the energy dispersion of exciton for different values of the total energy gap and 
interlayer separation and for the different dielectrics between the graphene layers. Results are presented 
for the parabolic approximation for the energy dispersion assuming the low-energy limit. When JC ~0.08 
nm -1 the parabolic approximation (|3U)) gives about 2% difference with respect to the exact numerical 
solution and this percentage decreases when JC deceases. The analysis of the results presented in Fig. 

shows that the energy dispersion decreases when the total energy gap increases. The same behavior 
can be observed for the dependence of the energy dispersion on interlayer separation: for given JC when 
interlayer separation increases the energy dispersion decreases (Fig. \Bp). In Fig. [__; is given the energy 
dispersion of exciton for the different dielectric placed between two graphene layers. As it seen, there is 
small increase for the energy dispersion: for a smaller value of the dielectric constant the energy dispersion 
of exciton becomes bigger. 



V. TWO BODY PROBLEM IN A GAPPED GRAPHENE LAYER 



Now we consider an electron and a hole located in a single gapped graphene layer with the energy gap 

parameter 8 and exciton formed by the electron and hole located in this graphene layer. Putting the gaps 

5i = 80 = 6 in ([5]). we obtain the Hamiltonian that describes two interacting particles in one graphene 

2 _____ 
layer interacting via Coulomb potential V(r) = 47T e eo ^ r ■ Using this potential we rewrite Eq. (|26|) in the 

form 



2e(hv F ) 2 V 2 r e 2 \ 
e 2 — 4<5 2 Answer J 



(Hv F JC) 2 



2c 

Eq. p5[) can be rewritten in the form of the two-dimensional hydrogen atom: 



(35) 



where 



e 2 \ 

-Ji(e)V 2 - cb aa = F {e)(j> aa , (36) 

4vT£ £r J 



2e(hv F ) 2 (hv F JC) 2 
1 = ' ^-^^le—- (37) 
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FIG. 2: Excitons in two graphene layers separated by the dielectric. The energy dispersion of exciton (a) for the 
fixed total energy gap and the different graphene interlayer separations, (b) for the fixed interlayer separation for 
different values of the total energy gap, and (c) for the fixed total energy gap, graphene interlayer separations 
and the different dielectrics between graphene layers. 
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The solution of the two-dimensional hydrogen atom equation (|30[) is well known [39h41| and is given by 



M*) = - 



4J 1 (e)e 2 (n-l/2) z 



(38) 



where n = 1, 2, 3, . . . are the quantum numbers, and the wave function in terms of associated Laguerre 
polynomials is given by 



<t>aa n i(r) = P 



(n- 1 



\l\ + n - (2n - 1) 



1/2 



(39) 



where /? = e 2 / [(n — l/2)eJ r i(e)], <p is the polar angle, L%(x) are the associated Laguerre polynomials, the 
quantum numbers I can take the values I = 0, ±1, ±2, . . . , ±(n — 1). 

After simplification Eq. (|38[) can be rewritten in the form of the following quadratic equation 

(C + 8 7 )e 2 - 4 7 fc 2 - ACS 2 = , (40) 

where 7 = (hvp) 2 , k = HvfK ~vpP, C = e 4 / (47re e(n — 1/2)) 2 . 
The solutions of Eq. (|40)) are given by 



\ 1/2 



e = 2 



-f(hv F K.y + CS 
C + 87 



(41) 
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Eq. (|4ip gives the energy-momentum dispersion e(/C) of the electron and hole that are bound via 
Coulomb potential in a single graphene layer. Since our interest is small kinetic energy, therefore, for 
small IC we expand Eq. (|4ip with respect to IC 2 : and approximate e(IC) by the first two terms of Taylor 
series 



E b 



(NCf 
2M 



where £5 is the exciton binding energy is given by 



E h = 26 



and M. is the effective mass of exciton given by 



C 



1/2 



(42) 



(43) 



(44) 



We note the exciton effective mass Ai increases when the gap 6 increases as it is shown in Fig. The 
result of calculation of the exciton energy dispersion in graphene layer for different gap energy 6 in low 
energy parabolic approximation is given in Fig. [3J). When IC ^0.08 nm" 1 the parabolic approximation 
(f4*2"j) gives less than 0.5% difference with respect to the exact solution (|4"Tj) and this value decreases when 
IC deceases. According to Fig. [Jt>, the exciton energy spectrum at the same quantum number n increases 
with the increase of the gap energy 5. However, for two graphene layers separated by the dielectric the 
energy dispersion decreases when the total energy gap increases. Also the comparison of the exciton 
energy distribution in graphene layer and in two graphene layers separated by the dielectric shows that 
the exciton energy in graphene layer always bigger than in two graphene layers separated by the dielectric. 




> 



Uj 



LU 



0.4 
0.35 

0.3 
: 0.25 
! 0.2 

2. 

5 0.15 
? 0.1 
0.05 




— 50 meV 

— 100 meV 

— 150 meV 

— 250 meV 







































































































































(a) 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 
K (nrrf 1 ) 

(b) 



FIG. 3: Excitons in the single graphene layer, (a) The effective exciton mass as a function of the energy gap.(b) 
The energy dispersion of exciton for the different values of the energy gap. 



VI. CONCLUSIONS 



Now let us return to the fundamental and practical question related to finding the solution of a problem 
of two interacting Dirac particles that form the exciton in gapped graphene layer and in two gapped 
graphene layers. In low-energy limit this problem can be solved analytically, and we obtained the energy 
dispersion and wave function of the exciton in gapped graphene layer and in two gapped graphene layers 
separated by a dielectric. The excitons were considered as a system of two Dirac particles interacting via 
a Coulomb potential V(r). In general case the center-of-mass and can not be separated from the relative 
motion even though the interaction depends only on the coordinate of the relative motion. This is caused 
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by the chiral nature of Dirac electron in graphene and the electron and hole locations and momenta are 
4-vectors and must satisfy mass-shell conditions. The analytical solution for the energy dispersion and 
wave function was obtained by introducing the transformation for the separation of the center-of-mass and 
relative motions for two particles in graphene that allows to reduce the Dirac-like equation for spinor to 
the Schrodinger-like second order differential equation for the component of the spinor. In the parabolic 
approximation for the energy dispersion found the effective mass of the exciton which are the function 
of the energy gap in the single graphene layer and function of the energy gaps and interlayer separation 
in the case of two graphene layers separated by the dielectric. Firstly, we can conclude that the exciton 
effective mass increases in the both cases for a single graphene layer and for two layers graphene as the 
gap energy increases. Also the exciton effective mass increases when the interlayer separation increases. 
Therefore, by tuning the energy gaps in graphene layers and changing interlayer separation one can 
get the desirable value for the effective exciton mass. This is very important for the system of many 
excitons when this system is considered as a dilute gas of excitons that forms Bosc-Einstcin condcncatc 
and undergoes to the Kosterlitz-Thoulcss phase transition to a supcrfluid phase. By decreasing the mass 
of the exciton one can increase the Kosterlitz-Thoulcss transition temperature. Secondly, for the exciton 
in graphene layer the energy dispersion increases with the increase of the gap energy. However, for the 
exciton in two graphene layers separated by the dielectric the energy dispersion decreases when the total 
energy gap increases as well as it decreases when the interlayer separation increases. 
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